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Abstract 

We present the results of a detailed calculation of the excitation spectrum of states 
with quantum numbers J^^ = O"*"^, 1 and 2"'""'" in the three-dimensional SU(2) 
Higgs model at two values of the scalar self-coupling and for fixed gauge coupling. 
We study the properties of Polyakov loop operators, which serve to test the con- 
fining properties of the model in the symmetric phase. At both values of the scalar 
coupling we obtain masses of bound states consisting entirely of gauge degrees 
of freedom (glueballs), which are very close to those obtained in the pure gauge 
theory. We conclude that the previously observed, approximate decoupling of the 
scalar and gauge sectors of the theory persists at large scalar couplings. We study 
the crossover region at large scalar coupling and present a scenario how the confin- 
ing properties of the model in the symmetric phase are lost inside the crossover by 
means of flux tube decay. We conclude that the underlying dynamics responsible 
for the observed dense spectrum of states in the Higgs region at large couplings 
must be different from that in the symmetric phase. 
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1 Introduction 



Over the past few years the SU(2) Higgs model in three dimensions has been the subject 
of many thorough numerical studies by means of lattice Monte Carlo simulations |]T| - |§] • 
The main incentive for these investigations was to clarify the nature of the electroweak 
phase transition in the framework of the dimensional reduction programme §], where it 
is much easier to obtain precise numerical results than in the full four-dimensional theory 
at finite temperature 0. In these studies it has been established that the electroweak 
phase transition is weakly first order for Higgs masses up to mn ~ 70 GeV, while for 
Higgs masses rriH > 80 GeV the transition disappears and turns into a smooth crossover 
P, 1^. Hence, the Higgs region and the confinement region in parameter space are indeed 
analytically connected, as was conjectured a long time ago for models with fixed-length 
Higgs fields 0. A review of the existing lattice results as well as a comparison with 



perturbation theory can be found in |10 



While the nature of the phase diagram and the order of the phase transition are 
now well determined, our understanding of the symmetric phase is still incomplete. Cal- 
culations of the mass spectrum based on the use of gauge-invariant operators lead to a 
picture of a confining symmetric phase with a dense spectrum of bound states. Recent 
attempts to describe the spectrum in terms of bound state models may be found in 



111 p!2| . This part of the phase diagram is of great intrinsic interest because it exhibits 



the same qualitative behaviour as QCD. In this region the non-Abelian gauge theory has 
confining properties, yet the presence of matter fields in the fundamental representation 
leads eventually to the screening of charges at large distances. This apparent similarity 
between QCD and the symmetric phase of scalar-gauge models opens the possibility to 
study some aspects of QCD in a model that is computationally much less demanding. 

In the Higgs phase, on the other hand, for small enough scalar couplings a pertur- 
bative calculation of the spectrum is applicable and leads to results in good agreement 
with those from lattice simulations. Due to the analytic connectedness of the phase di- 
agram all mass eigenstates in the Higgs region may be continuously mapped into their 
counterparts in the confinement region. This offers the possibility to study the onset of 
non-perturbative physics and confinement by smoothly moving from the Higgs region 
into the confinement region of the phase diagram. 

In this paper we elaborate on our previous computation of the mass spectrum of 
the SU(2) Higgs model in 2+1 dimensions]^]. Here, however, we move away from the 
original context of the electroweak phase transition, and our main interest lies in the 
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dynamics of the three-dimensional theory and what it might have in common with QCD. 
In particular, we want to obtain a more complete picture of the excitation spectrum and 
investigate in more detail the previously observed approximate decoupling of the scalar 
and gauge sectors of the model 0]. Thus, in addition to the 0+"*" and 1 states considered 
before, we now also compute correlations of 2"^^ and Polyakov loop operators. Although 
our analysis does not yet incorporate multiparticle states, it leads to a more complete 
picture of the symmetric phase. 

In contrast to our previous work, where we carefully studied the approach to the 
continuum limit, we now consider the theory at a fixed value of the gauge coupling, 
i.e. (3g = 9, where our earlier results show that the physics is already very close to 
the continuum limit. Furthermore, at large scalar coupling it has been demonstrated 
numerically [@] that the disappearance of the phase transition survives in the continuum 
limit and that no significant lattice artefacts are to be expected for jSc = 9. 

Our main results are as follows. The previously observed approximate decoupling of 
states composed purely out of gauge degrees of freedom in the 0"*"+ channel extends to 
the 2~^~^ channel. The masses of these states are very close to their values in the pure gauge 
theory. Furthermore, the decoupling occurs at both small and large values of the scalar 
self-coupling, and in the latter case it can be observed also deep in the crossover region. 
Second, at large scalar coupling we find that a dense spectrum of states is observed even 
in the Higgs phase, yet the underlying dynamics must be quite different, as confinement 
is no longer observed. We present a scenario how the confining properties of the theory 
are lost inside the crossover region: within the QCD context it is natural to use the 
concept of fiux tubes in the symmetric phase as a manifestation of confinement. As 
one enters the crossover region the fiux tube develops a rapidly growing decay width and 
eventually becomes so unstable that there is no longer any remnant of linear confinement 
(i.e. a range of distances where the static potential shows a linear behaviour). Finally, 
the inclusion of correlations of Polyakov loops changes some of our previous conclusions 
about the infiuence of finite-volume effects. In particular, we find that finite-size effects, 
especially in the lowest 0"''^ state are less severe than was claimed before H. 



Some of the results presented in this paper have been summarised elsewhere [13, 14 



We also draw the reader's attention to some recent related work 15 



The remainder of this paper is organised as follows. In section 2 we briefly summarise 
the lattice techniques used in this work. Section 3 contains a discussion and comparison 
of the properties of the Higgs and confinement phases for small scalar self-coupling. In 
section 4 we describe how the properties of the Higgs and confinement regions of the 
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phase diagram change when the size of the scalar self-couphng is increased. In sections 
we follow the analytic connection of the mass spectrum through the crossover. Results 
indicating flux tube decay in the crossover region are presented in section 6, and section 7 
contains our conclusions. 



2 Methodology 

We work with the three-dimensional lattice action given by 

S[U,<P] = /?GE(l-^Trf/,)+5:j-/?^,X^iTr(0t(a;)f/^(x)0(x + /i) 

p ^ ^ ^ X [ fj.=l ^ 

+ iTr(0t(a;)0(a;)) + (3n [^Tr (0t(x)0(a;)) - iH, 



(1) 



where U^{x) G SU(2) is the link variable, Up denotes the plaquette, and is the 
scalar field. The bare parameters Pg, Ph and Pr are the inverse gauge coupling, the 
scalar hopping parameter, and scalar self-coupling, respectively. In our simulations we 
fix the combination of parameters 

9i~ ' 

where Xs/g^ is the ratio of the scalar self-coupling and the gauge coupling in the con- 
tinuum formulation. Exact relations between the lattice and the continuum parameters 
at the two-loop level governing the approach to the continuum limit can be found in 

i) 



161 . In this work we consider the cases of small, i.e. A3 <^ and large, i.e. A3 ~ 0{g 



scalar self-coupling, by choosing X^/g^ = 0.0239 and Xs/g^ = 0.2743, respectively. At the 
small value, which was used in our earlier study [Q, the confinement and Higgs phases 
are separated by a first order phase transition, whereas for the larger value of X^/g^ the 
transition has turned into a smooth crossover^. 



2.1 The blocking procedure 

In any mass calculation on the lattice it is important to have a good projection property 

of the interpolating operators in order to obtain a reliable signal. For gauge theories 

^Note that at tree level, the ratio As/yf is proportional to the ratio of Higgs and M^-boson masses in 
the four-dimensional theory at finite temperature. As an orientation for those readers who are familiar 
with studies of the electrowcak phase transition, we add that in the framework of dimensional reduction, 
our values of Xs/g^ correspond to tree- level, zero temperature Higgs masses of 35GeV and 120 GeV, 
respectively [|l[ . 
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with or without scalar fields it has been demonstrated that the projection of the basic 
operators onto the desired state can be considerably enhanced by employing so-called 



"blocking" or "fuzzing" techniques [17, 18, 19 



As in our previous work[Q, we construct blocked link variables U^\x) at blocking 
level n according to the procedure originally defined in[|r^. Furthermore, a blocked 
version of the scalar field, at blocking level n is constructed according to 




(3) 

Note that we will take correlations in the 3 direction and so our blocking always remains 
within the (1, 2)-plane. In addition, we consider blocked scalar fields (f)^"''^\x) with non- 
local contributions from one spatial direction j only: 
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+ f/f-i)t(a;-j)0(«-i.^-)(x-j)}, J = 1,2. (4) 



It is obvious that all blocking procedures employed are gauge invariant and can be 
iterated in order to create a large set of operators of various spatial extensions in an 
efficient manner. In practice, the blocking levels have to be tuned in order to identify the 
operators with the best projection properties. 



2.2 Basic operators 

We now proceed to list the basic operators used in our simulation. Since we use a larger 
set of operators compared to our previous calculation, we have decided to change our 
notation slightly for clarity. Although in general we use blocked links and scalar fields, 
we suppress the superscripts which label the blocking level in this subsection, except 
when stated explicitly. 

After decomposing the scalar field (j){x) as 

(f){x) = p{x)a{x), p{x) > 0, a{x) G SU(2), (5) 

we consider the following basic operators involving only scalar fields or combinations of 
scalar and gauge degrees of freedom: 

Rix) = iTr(0t(a;)0(x)), 
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Lj{x) = -Tr{a\x)Uj{x)a{x + j)), j = 1,2 

V;{x) = ^Tr(rV(x)C/,-(x)a(x+j)), j = 1, 2, (6) 

where is a Pauli matrix. Note that at this stage we have not yet specified the quantum 
numbers J^*" of the operators. This is postponed to the end of this subsection when all 
basic operators have already been defined. 

We also consider Wilson loops Cij, i.e. gauge invariant operators constructed out of 
link variables only, e.g. 

Cl;\x) = Tt{u,{x)Uj{x + i)uUx + i)U]i^)) , hj = 1,2, ^ 7^ J, (7) 

and in addition to the elementary plaquctte C^^^, we also consider squares of size 2x2 
as well as rectangles of size 1x2, 1x3, 2x3. Hence, we have five versions of Wilson 
loop operators, viz. 

/^Ixl /~i2y.2 /^lx2 /^lx3 /^2x3 (q\ 

As we shall see, a larger set of this type of operator will significantly improve our ability 
to extract information about the gauge sector of the theory. 

An operator which is useful to probe the confining properties of a theory is the 
Polyakov loop, i.e. a loop of length L of link variables that winds around the spatial 
dimensions of the lattice, viz. 

Pf\x) = Tr n + ^J), 3 = 1, 2. (9) 

m=0 

The expectation value of the Polyakov loop vanishes if the theory is confining, and it 
projects onto the state consisting of a chromoelectric flux loop that closes through the 
boundary in the j direction. Correlations of P^^^ can then be used to extract a string 
tension from the lowest mode of its exponential fall-off, the so-called "torclon" mass. 
However, in a theory with matter fields, such as OTirs, one eventually expects the string 
to break beyond some large separation, due to pair creation (just as in QCD). In our 
model we also have a phase with spontaneously broken symmetry, the Higgs phase. Here, 
this interpretation of the Polyakov loop correlator no longer holds, the Polyakov loop 
has a large expectation value, and the torclon mass cannot be used to extract a string 
tension. Instead, since a weak coupling expansion of the Polyakov loop correlator applies 
in the Higgs phase, its leading term corresponds to an exchange of two spin-1 particles 
(i.e. H^-bosons), and therefore the lightest mass above the vacuum should be interpreted 
as a two-PF scattering state. 
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Apart from being relevant for the study of confinement, the Polyakov loop also plays 
an important role in understanding apparent finite-size effects in glueball calculations in 
the pure gauge theory . Since we are interested in the detailed structure of the spec- 
trum, and, in particular, higher excitations, we need to address the influence of torelons 
in the various spin channels. Although single torelons, arising from correlations of the 
Polyakov loop Pj^^ cannot contribute directly to correlations of operators C^^^, . . . , C^^^, 
they may well be relevant for the operators R{x) and Lj{x). This is because of the it- 
erative, non-local nature of the blocking procedure, such that highly blocked versions 
of R{x) and Lj{x) may contain contributions which, for periodic boundary conditions, 
wind around the lattice and may thus have a sizeable projection onto Pj^^ ■ Furthermore, 
torelon-antitorelon pairs are known to give rise to large finite-volume effects in glueball 
calculations |2^. Apart from studying correlations of Polyakov loops per se, we have de- 



cided to construct single-torelon and torelon-pair operators in the 0"'""'" and 2"*"+ channels 
and to study their correlations as a safeguard against finite-size effects. 

We now give a complete list of the types of operators with quantum numbers J^*" 
and also introduce the symbols which will later label the contributions from the various 
operators to a given state. 

0++ channel: 
R: R 
L: Li + L2 

C: symmetric combinations of C^^\ C^^^, C^^^, C^^^, C^^^ 
P: P[^^ + P^^^ 
P.: Pi'' ■ 



T: {pi'^y + {pi'^y 



1 channel: 

V: V;, J = 1,2 a = 1,2,3 

2^+ channel: 

R: = i {Tr[0("'i)t0(".i)] _ Tr[0(".2)t0(n,2)]| ^ ^ > ^ 

Li — L2 

antisymmetric combinations of C^^^, C^^^, C^^^ 



L 
C 
P 
T 



p{L) _ p{L) 



p(.))^_(p(.)y 
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Here, the operators are understood to be zero-momentum sums of the corresponding 
operators defined in eqs. The labels R, L refer to operators containing scalar 

fields, whereas the Wilson loop operators C are composed only out of gauge degrees of 
freedom. The labels P, Pd denote torelon operators, and T labels pairs of torelons in the 
respective channels. Note that the definition of the operator b!"^^ in the 2^^^ channel is 
only possible with the kind of blocking defined in eq. (^). 

Our list of operators does not contain multiparticle operators, and therefore we do 
not have control over continuum states in the spectrum^ Here we wish to remark that we 
have evidence that continuum states have only small overlaps on the basis of operators 
used in this work, so that we are still able to extract the bound state mass spectrum 
reliably. 

2.3 Matrix correlators 

In order to compute the excitation spectrum of states with given quantum numbers, we 
construct matrix correlators by measuring all cross correlations between different types of 
operators at several blocking levels @] . This correlation matrix can then be diagonalised 
numerically following a variational method. For a given set of N operators 0i we find the 
linear combination that minimises the energy, corresponding to the lightest state. 

The first excitation can be found by applying the same procedure to the subspace 
which is orthogonal to the ground state. This may be continued to higher states so 
that we end up with a set of eigenstates $j, z = 1, . . . , given by 

N 

^i = Yl ^ik(t)k- (10) 
k=l 

The coefficients quantify the overlap of each individual operator </)fc used in the simula- 
tion onto a particular approximate mass eigenstate $j. For a complete basis of operators 
this procedure is exact. In practice, the quality of the approximation clearly depends 
on the number of original operators and their projection properties. The determined 
eigenstates are removed from the basis for the higher excitations, so that the basis for 
the latter gets smaller and the corresponding higher states are determined less reliably. 

For our main spectrum calculation we work with iV = 24 operators in the 
channel, i.e. (pk ^ R-, L, C, P, Pd, T, each considered for at least two blocking levels. In 
the 2"*"+ channel we choose = 23 and (pk € R, L, C, P, T. In the symmetric phase this 
•^Notc that on a finite volume, those states have a discrete spectrum. 
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typically enables us to obtain mass estimates for the first 10-12 states in each of these 
channels. 

For spin-1 particles, we construct a 4 x 4 matrix correlator from the operator of 
type V considered at four blocking levels. Thereby we can extract masses for the first 
two or three states only. In the computation of separate correlations of the Polyakov loop 
operator itself, we also compute a 4 x 4 matrix correlator using different blocking levels. 

For a more detailed discussion of the variational method based on matrix correlators, 
along with numerical investigations to test its efficiency, we refer to our previous work @] . 



2.4 Simulation and analysis details 

The Monte Carlo simulation of the lattice action in eq. (|I]) is performed using the same 
algorithms as in our previous work. The bare parameters Pg, Ph and are fixed by 
applying the constraints Xs/g^ = 0.0239 and Xs/g'^ = 0.2743, respectively, and using 
the same procedure as in section 3 of ref . . For the update of gauge variables we use a 



combination of the standard heatbath and over-relaxation algorithms for SU(2) ^ 
The scalar degrees of freedom are updated using the algorithm described in[^. In 
general, the applicability of this algorithm is restricted to cases where the self-coupling 
Pr is not too large, as otherwise the acceptance rate decreases substantially. In our 
simulation, however, we did not observe any significant drop in the acceptance rate, even 
at our larger value of Pr. Thus the algorithm may be applied safely in the region of the 
parameter space explored in this paper. 

As before, we define a "compound" sweep to consist of a combination of one heatbath 
and several over-relaxation updates of the gauge and scalar fields. In the Higgs and 
confinement regions of the phase diagram at both small and large scalar coupling we 
typically gathered about 12000 compound sweeps for our main spectrum calculation. In 
the crossover region we also monitored the spectrum at several, closely spaced values of 
Ph, for which about 5000 compound sweeps were accumulated each. Thus, the statistics 
in our present study is much lower than in our previous work. At first sight this may seem 
surprising, since we now seek to obtain a more detailed structure of the mass spectrum. 
However, we wish to remark that our basis of operators is much larger than in , so that 
a loss in statistics is compensated for by using the information from a greater number of 
correlators in a particular channel. In addition, the diagonalisation procedure requires a 
large basis of operators in order to yield reliable results for the excitation spectrum. 

Our standard lattice sizes in this study are 36^ and 24^, which, on the basis of our 
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previous study, do not lead to significant finite-size effects in the 0"''^ and 1 channels. 
In order to study finite-volume effects in more detail, especially in the 2"'""'" channel, we 
have also simulated lattices of size 16^ ■ 24. 

Our mass estimates are extracted from two-parameter, correlated, single-exponential 
fits over a finite interval [ti, ^2] to correlation functions of eigenstates $j (c.f. subsec- 
tion 3.4 in ref. Q). We have checked explicitly that the alternative fitting formula moti- 
vated in p4|, where an additional constant is considered, does not produce an appreciable 



difference in our mass estimates. 

In cases where the effective masses calculated from the correlation function did not 
show a plateau long enough to perform a correlated fit, we used effective masses for 
our best estimates. This was also our preferred method for the scan of the crossover 
region, where we were more interested in qualitative features rather than very precise 
mass estimates. Results which were obtained from effective masses, or from eigenstates 
whose individual overlaps aik were small throughout, are marked by an asterisk in our 
data tables in sections^ and^. 

Statistical errors are estimated using a jackknife procedure for which the individual 
measurements have been accumulated in bins of 500 sweeps. For some of our most precise 
results we also quote a systematic error, which was obtained by quoting the difference 
between our best estimate and the result from an uncorrelated fit, and, in some cases, 
the result obtained using an alternative fitting interval. 



3 Results at small scalar coupling 

To explore small scalar coupling physics we fix X^/g^ = 0.0239 as in ref. For this 
value the system exhibits a strong first order transition upon variation of jSn and the 
Higgs and confinement properties on each side of the transition are very pronounced and 
easily distinguishable. In the symmetric phase we choose the point Ph = 0.3438 and in 
the Higgs phase = 0.3450. The O"*"^ and 1 channels of the mass spectrum at these 
parameter values have already been determined in our previous work . Here we repeat 
this calculation with a larger basis of operators in order to obtain more details about the 
excitation spectrum. We do so at (3g = 9, where we know from our previous study that 
we are already close to the continuum limit of the theory. 
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Figure 1: (a) The lowest states of the spectrum in the confinement (left) and Higgs (right) 
region at small scalar coupling, /5g = 9. The lattice sizes are 36^ in the confinement and 
24^ in the Higgs phase, respectively. Full symbols indicate states that receive predomi- 
nantly contributions from Wilson loops, (b) The coefficients aik of the 24 operators used 
in the simulation on 36^ for the five lowest 0^"^ eigenstates in the confinement phase. 

3.1 Mass eigenstates in the Higgs and confinement regions 

The lowest states of the mass spectrum at the two points investigated are shown in 
Fig.[^(a). In the confinement region one observes the now famihar, dense spectrum of 
states for the three spin channels. As an illustration of the contributions of different 



types of basic operators to the mass eigenstates, we plot the coefficients an. (c.f. eq. ([Tol)) 
in Fig.|l|(b) for the lowest five eigenstates in the 0"^^ channel in the confining phase. 

Our main results in the Higgs phase are listed in Table|l|. Here, the only low-lying 
states are the familiar Higgs and ly-boson states, which are also well described by per- 
turbation theory. Note in particular that the lowest mass in the 2"*"+ channel is rather 
large and consistent with that of a VT-scattering state - as are the excited states in 
the 0"'""'" channel. In particular, consider one of the first excited states in the 0"^"*" channel. 
This state has almost exclusively contributions from Wilson loops, i.e. it consists mainly 
of gauge degrees of freedom. This suggests interpreting it as a ly-scattering state, and 
Table ^ shows that the mass of this state is indeed compatible with twice the ly-mass. 
This result is also expected from perturbation theory, where the term with the slowest 
exponential fall-off corresponds to a two-W^ exchange diagram. It is natural to ask at this 
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Table 1: Results for all three spin channels in the Higgs region at (3h = 0.3450. The sec- 
ond error in some results is an estimate of systematic errors as described in subsect. 

point why we do not see a two-Higgs state as the first excitation, since such a state is 
hghter than the two-W^ state. The reason is that our operators appear to have practically 
no projection on such a state. In order to see a two-Higgs state it would be necessary to 
study the correlations of an operator which by construction couples to such a state in 
leading order, for instance {R{p = 0) — (-R))^- This, however, is beyond the scope of this 
work. 

Our main results in the symmetric phase are listed in Tables ^|-|^. Here and in the 
following tables, the ordering of states $j was obtained during the diagonalisation proce- 
dure from the effective mass of the corresponding correlator on the first timeslice. With 
our larger basis of operators we are able to identify a large number of states (e.g. up to 16 
states in the O"*"^ channel), and quite a number of them receive their dominant contribu- 
tions from Polyakov-loop or torelon-pair operators. Note in particular the change in the 
operator content in the ground state of the 2"*"+ channel as the lattice size is increased. 
This will be discussed in more detail in connection with finite-size effects in the next 
subsection. Our results for the lowest masses are in agreement with our earlier study 
(c.f. Tables 2 and 3 in |^), although we now obtain a slightly lower value for the mass in 
the 1 channel. 

The physical states consist of a dense spectrum of bound states. The second excited 
state in the 0"^"*"-channel is a pure gauge excitation which, in this phase, cannot be per- 
turbatively related to a scattering state of two spin-1 particles. Instead it was suggested 
in 1^ to interpret it as a glueball in analogy to the glueballs of pure gauge theory. This 
interpretation is based on the observation that the state is composed predominantly 
from gauge degrees of freedom and that its mass nearly equals that of the lightest scalar 
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Table 2: Mass spectrum and dominant contributions from the operator basis in the 0~^^ 
channel at small scalar self-coupling in the symmetric phase. 
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r 1 

G 


$8 
$9 


1.47* 
1.49(4) 


c 

R 






16^ • 24, 




= 0.3438 




No. 


aM[2-+] 


Ops. 


No. 


aM[2+-] 


Ops. 


$1 


0.361(12) 


P 


$6 


1.5(1)* 


C 


$2 

$3 


0.843(27) 
0.86(2) 


T 
L 


$7 

$8 


1.4(1)* 
1.55(4) 


L,R 
C,T 


$1 


0.98(4) +5 


P 


$9 


1.6(2)* 


L 


$5 


1.18(5) tl 


C 


$10 


1.79(7) 


C 



Table 3: Mass spectrum and dominant contributions from the operator basis in the 2++ 
channel at small scalar self-coupling in the symmetric phase. 



pH = 0.3438 




L'^ .T = 36^ 


L'^ ■T = 24^ 


• T = 16^ ■ 24 


No. 


aM[l—] 


aM[l—] 


aM[l--] 


$1 


0.595(5) 


0.582(14) 


0.585(15) 


$2 


0.857(8) 


0.863(9) +° 


0.99(2)* 


$3 


0.98(8)* 


1.068(18) 


1.15(4) 



Table 4: Mass spectrum in the 1 channel at small scalar self-coupling in the symmetric 
phase. 



13 



glueball in the pure SU(2) gauge theory at this value of [3g- As Table^ shows, there is a 
corresponding state in the 2^"'" channeL We also find higher excitations of these lowest 
glueball states for both quantum numbers. 

In Table ^ we compare some properties of the Polyakov loop operator in the Higgs 
and confinement regions. As explained in subsection|2j^, the Polyakov loop has a large 
expectation value in the Higgs phase, indicating the absence of confinement. Similar to 
the correlator of Wilson loops, the perturbative expansion of the Polyakov loop correlator 
for weak coupling is dominated by a iwo-W exchange diagram, and indeed we find its 
exponential fall-off to be compatible with twice the W-m.a,ss. These findings are supported 
by the observation that the mass extracted from the correlation of Polyakov loops is 
roughly independent of the spatial lattice sizeL, in sharp contrast to situations where 
confinement is observed. 

In the symmetric phase, on the other hand, we find that the expectation value of 
the Polyakov loop is consistent with zero, indicating confinement. Here, ordinary weak 
coupling expansions fail, and the exponential fall-off of the correlation of Polyakov loops 
is incompatible with twice the mass of the spin-1 state, c.f. Table |^. Instead, one interprets 
it as the mass of a loop of chromoelectric flux, aMp, which can be related to the string 
tension ai 

{P^''\x)pj''^\0)) ^ e-''*^^^^)*, aMpiL) = a^aLL. (11) 



Here L is the spatial length of the lattice and hence of the flux loop. An estimate for the 
string tension in inflnite volume is then provided by the relation p5| 

aVoo = aVi + (12) 

We add that our numerical values for the string tension a^fa^ shown in Table § are only 
around 3% smaller than in the pure gauge theory | 



3.2 Finite-volume effects and toroidal operators 

In this subsection we consider the issue of flnite-size effects in more detail. In particular, 
we want to clarify the role of toroidal (or winding) operators, i.e. operators constructed 
out of Polyakov loops, for flnite-volume studies in QCD-like situations. This is of par- 
ticular relevance in the symmetric phase of our model, where we have reported strong 
flnite-size effects in the ground state of the O^"*" channel 0]. Although we include this dis- 
cussion in the section on small scalar coupling, our flndings are relevant at large coupling 
as well. 
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Confinement 
I3h = 0.3438 


Higgs 
/3h = 0.3450 


L'-'-T 


36^ 


243 


16^ ■ 24 


16^ ■ 24 


243 


{Pi") 

aMp 

aM[l—] 


0.001(10) 
0.879(6) +° 
0.595(5) +° 
0.1575(5) +° 


0.001(10) 
0.570(4) 

0.582(14) 

0.1570(6) +° 


0.14(2) 
0.366(3) 
0.585(15) 
0.1579(5) 


6.535(6) 
1.67(4) 

0.839(4) 


6.724(5) 
1.54(4)* 

0.837(3) +° 



Table 5: Properties of the Polyakov loop operator in the confinement and Higgs region 
for = 9. (Pi^^) here refers to the unsmeared operator. 



Toroidal operators, such as the Polyakov loop, project mainly onto states encircling 
the periodic boundary conditions, sometimes called "torelons" |2^ (see also eq. ([Til)). 
These torelon states can interfere with the mass spectrum of the theory in infinite volume 
if 

aVLL<aM, (13) 

where M is the mass of an eigenstate of the Hamiltonian. It should be stressed at this 
point that torelon states exist in the continuum limit for finite box size and periodic 
boundary conditions. Thus they are not lattice artefacts, i.e. artefacts arising from the 
discretisation of the theory. 



As explained in subsection |2.2| , we use blocked operators in the 0"^"^ and 2++ channels, 
which, due to the iterative nature of the blocking procedure, may contain contributions 
which wind around the lattice and can thus have a sizeable projection onto torelon states. 
If relation ([T3| ) applies this can then lead to a possible misidentification of the ground 
state and/or higher excitations of the Hamiltonian. This scenario does not apply in the 
1 channel, where it is impossible to construct a torelon operator with the same quan- 
tum numbers. Furthermore, the overlap between single torelons and glueball operators 
will be completely suppressed in the pure gauge theory where we have confinement and 
the associated Z2 symmetry, under which torelon and glueball operators transform dif- 
ferently. However, torelon-antitorelon pairs can have a big influence on the spectrum of 
glueballs and are therefore included in our study. 

As an illustration of the issue, consider Fig. |1| (b) and Table 0. The analysis of the 
operator content shows $4 to have almost exclusively P-content. The mass associated 
with it is a flux loop mass appropriate to extract the string tension, but it is not part of 
the spectrum of the Hamiltonian in infinite volume and has to be removed. An example 
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how these toroidal states mix with the physical ones and have to be carefully disentan- 
gled is given for the case of the 2"'"+ spectrum in Fig. ^ and Table ^ On the smaller 




Figure 2: The coefficients of the 23 operators used in the simulation for the five 
lowest 2++ eigenstates, (a) L = 24, (b) L = 36. 

lattice, Fig. |(a), the lowest state is a mixture of R/L and P-type operators, with a 
slight dominance by P. Quoting the corresponding mass as a physical state close to its 
infinite-volume limit would be erroneous, as is elucidated by considering the same opera- 
tors on a larger volume. The masses arising from winding operators have become heavier 



(see eq. (|T3|)), while the non-winding operators produce mass values which hardly change 
under variation of the volume, and are thus already close to infinite- volume physics. Com- 
bining the information about the overlaps of individual operators with their behaviour 
as the spatial lattice size is increased makes the identification of the masses which are 
part of the spectrum in the infinite-volume limit unambiguous. 

As a result, some of the conclusions concerning the finite- volume behaviour in ref. 0] 
have to be revised. First, the large finite- volume effects in the O"*""*" channel clearly arise 
from toroidal operators, as the relation ([T3|) certainly holds in the 0"'""'" channel for the 
smallest lattices considered in our previous work. Second, our new findings remove the 
anomalous behaviour of the first excited state. This can be understood since the ordering 
of the excitation spectrum changes with increasing lattice size, as toroidal operators will 
produce more and more massive states. 
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-Higgs 


Pure Gauge 


Gauge- 


-Higgs 


Pure Gauge 


Scalar 


glueball 


glueball 


Scalar 


glueball 


glueball 


u.oCdyoj \j 






U. ( 04(^ i ) 






0.706(5) +° 






0.99(3) 






0.95(3) 


0.74(3)* 


0.767(6) 


1.28(2) 


1.23(2) 


1.26(2) 


0.92(3) tl 






1.25(15)* 








1.08(4) tl 


1.08(2) 




1.46(2) 


1.50(5) 




1.35(2) 


1.27(2) 


1.67(6) 


1.64(4)* 


1.77(6) 



Table 6: Final mass estimates using the data obtained on 36^ at Pg = Ph = 0.3438 in 
the symmetric phase at small scalar self- coupling. Our data for the glueball are compared 
to results obtained in the pure gauge theory. 

3.3 The physical spectrum 

After removing finite size effects introduced by considering toroidal operators, we are 
ready to present our final mass estimates for the spectrum in Table ^ In addition to the 
mass spectrum extracted from our simulations of the SU(2) Higgs model, the glueball 
spectrum of the pure SU(2) gauge theory is also given p7[ . 

For all cases considered we find that the glueball masses deviate from those in the 
pure gauge theory at the percent level at most. Thus, the glueball spectrum of the pure 
gauge theory is almost identically repeated in the scalar-gauge theory and appears to 
be entirely insensitive to the presence of the scalar fields. In addition, and apparently 
disjoint from this part of the spectrum, there are bound states of scalars which also have 
some gauge content. Combining this observation with the practically identical value of 
the string tension, we confirm our earlier conclusion that the pure gauge sector of the 
theory decouples almost completely from the scalar sector at this point in the symmetric 
phase. 

4 Results at large scalar coupling 

It is now interesting to study whether a significant increase of the scalar self coupling 
destroys the decoupling of the pure gauge sector observed for small couphng. For this 
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Figure 3: The lowest states of the spectrum in the confinement (left) and Higgs (right) 
region. Full symbols indicate states that receive predominantly contributions from Wilson 
loops. 

purpose we fix the ratio of scalar and gauge couplings to be about ten times larger 
than previously, Aa/^fg = 0.2743. At this parameter value the first order phase transition 
between the Higgs and the confinement regions upon variation of jSn has disappeared 
and turned into a smooth crossover p, |^. One would expect that the absence of a phase 
transition weakens the specific Higgs and confinement properties on the respective "sides" 
of the crossover. For the simulations the points jSn = 0.35262 and jSn = 0.35542 are 
chosen to represent the confinement and Higgs sides of the phase diagram, respectively. 

The lowest states of the mass spectrum are shown in Fig. |^. Higher excitations as 
well as the composition of all states can be found in Tables |7|-|TI|. The spectrum in the 



confinement region looks qualitatively the same as in the small scalar coupling case. 
In particular, we observe the same occurrence of glueball states in the O^"*" and 2"*"+ 
channels as well as their non-mixing with the bound states of scalars as for small scalar 
coupling, c.f. Tables ^, ^. By comparing the numerical values of the glueball masses at 
small and large scalar coupling, one finds that all glueball masses agree quantitatively, 
i.e. are statistically compatible, with those at small scalar coupling. Thus the previously 
stated decoupling of the pure gauge sector from the Higgs part of the theory extends 
over an order of magnitude increase in the ratio of scalar coupling to gauge coupling! 

On the Higgs side of the crossover, a dramatic change in the nature of the spectrum 
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L2.r = 


363, 


Ph = 


0.3526 




No. 


aM[0++] 


Ops. 


No. 


aM[0++] 


Ops. 


$1 
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R,L 


$8 
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S 


$2 


0.645(5) +° 


R, L 


$9 


1.33(3)(1) 


c 


$3 
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C 


$10 


1.21(9)* 


L 


$4 


0.850(7) +1 


P 


$11 


1.40(3)* 


P 
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0.785(19) 


L 




1.56(4)* 
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$6 


0.92(3) +° 


R,L 
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1.60(4) 


C 
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1.08(3) 1° 
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0.3526 




No. 
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No. 
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Ops. 


$1 
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0.324(4) 
0.543(9) 
0.661(7) 


L 
P 

R,L,P 


$9 

$10 
$11 


1.10(3) 
1.17(3) 
1.18(3) 


C,T 
C,T 


$4 
$5 


0.737(12) 
0.752(23) 


C 
L 


$12 
$13 


1.42(3)* 
1.3(2)* 


T 
L 


$6 


0.968(12) 


R,L 


$14 


1.53(3) 


C,T 


$7 
$8 


0.98(3) 
1.03(2) 


C,Pd 
C,Pd 


$15 


1.53(4) 





Table 7: Mass spectrum and dominant contributions from the operator basis in the 0'^'^ 
channel at large scalar self-coupling in the symmetric phase. 
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Ops. 
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Ops. 


$1 


0.568(5) tl P 


$6 


1.29(2) 


C,T 




0.734(23) 


L 


$7 


1.22(2) +0 


L 


$3 


1.02(2) 


R 


$8 


1.48(4) +0 


C 


$4 
$5 


1.12(2) 
1.22(3) 


P 

C,T 


$9 
$10 


1.40(15)* 
1.65(5) 


R 
R 



Table 8: Mass spectrum and dominant contributions from the operator basis in the 2 
channel at large scalar self-coupling in the symmetric phase. 



has taken place. The specific "Higgs-like" feature of having only the Higgs and VT-particle 
as low lying states with a large gap to higher excitations has entirely disappeared. Instead, 
we now have a dense spectrum of states on the Higgs side as well, which differs only 
quantitatively from the one on the confinement side. Note that the value for the scalar 



coupling considered here is so large that perturbation theory is not reliable anymore [28 



Let us now look at the expectation value of the Polyakov loop. The points where we 
analysed the mass spectrum correspond to the endpoints in Fig. ^, where the expecta- 
tion value of the Polyakov loop is shown. Note that now there is a small but non-zero, 
measurable value for the expectation value of the Polyakov loop on the confinement side. 
This is a first indication that colour charges are screened at large distances. On the other 
hand, at the Higgs end we measure a large expectation value, just as in the case of small 
scalar coupling. Comparing the two sides in Fig. ^ one might be inclined to conclude 
that the similarity of the spectra indicates an extension of the confining physics over to 
the Higgs side of the crossover, which seems not unreasonable given the absence of a real 
transition. However, from the large difference in the expectation value for the Polyakov 
loop we are led to conclude that the dynamics responsible for the dense spectrum on the 
Higgs side must be very different from that on the confinement side. In particular, no 
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Table 9: Mass spectrum and dominant contributions from the operator basis in the O"*""*" 
channel at large scalar self-coupling in the broken phase. Labels in brackets indicate that 
the contribution from the corresponding operator are small but significant. 
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Table 10: Mass spectrum and dominant contributions from the operator basis in the 2++ 
channel at large scalar self-coupling in the broken phase. Labels in brackets indicate that 
the contribution from the corresponding operator are small but significant. 
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Table 11: Mass spectrum in the 1 channel at large scalar self-coupling in the symmetric 
and broken phases. 
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Table 12: Final mass estimates using the data obtained on 36^ at = 9, (3h = 0.3526 in 
the symmetric phase at large scalar self- coupling. Our data for the glueball are compared 
to results obtained in the pure gauge theory. 

flux loops exist on this side of the crossover. This may be an indication that what we 
see in this region of the phase diagram is a strong scalar coupling Higgs regime. More 
numerical and analytic evidence is required, though, to verify this interpretation. 

Our final mass estimates for the physical states in the confinement phase, with the 
toroidal operators removed as in the previous section, are given in Table |l^. 

5 The crossover region 

As we have seen in the last section, the mass spectra of the Higgs and confinement re- 
gions are rather similar at large scalar coupling, although the underlying dynamics must 
be quite different. Since there is no phase transition but only a smooth crossover "sep- 
arating" the respective regions at large scalar coupling one can study how the spectra 
are continuously connected. Apart from mapping the spectra in the Higgs and confine- 
ment regions onto each other, this also offers the opportunity to follow the change in the 
dynamics as one moves from one regime into the other. 

We start by considering the expectation values of the Polyakov loop and the length of 
the Higgs field, displayed in Fig. |. In most simulations dealing with the phase transition 
R has been used as an "order parameter" , as it jumps across a first order phase transition 
from small values in the confinement phase to large values in the Higgs phase. As Fig. ^ 
illustrates, this jump has entirely disappeared, the change in (R) is smooth and very 
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Figure 4: (a) Vacuum expectation values for R and P. (b) Susceptibility of a smeared 
R- operator. 

slow, with its values differing by only ~ 30% at the end points studied in the last section. 
On the other hand, as mentioned in the previous section, the expectation value of the 
Polyakov line is quite different at both ends of the curve. Because it probes confining 
behaviour and the change in dynamics in this region of parameter space, the expectation 
value of the Polyakov loop is more sensitive to the changes between the regions and 
hence makes a better "order parameter". However, strictly speaking it has to remain 
non-zero on the confinement side as well, so the change is smooth and quantitative only, 
as expected for a crossover. 

Despite the smoothness of the crossover, it is possible to define a critical hopping 
parameter through the peak of the susceptibility of various quantities. We obtained the 
best signal for the susceptibility of a smeared version of our i?-operator, 

X = {R') - {Rf . (14) 

In Fig. ^ (b) one easily identifies a pronounced peak in the susceptibility indicating large 
fluctuations in the operator under consideration. As in ref. this peak serves to locate 
the critical value of the hopping parameter f3fj at which the crossover takes place, and 
thus still represents a line "separating" the Higgs and confinement regions even in the 
absence of a phase transition. At /?g = 9 we estimate it to be P'^ = 0.35406(5). 

The variation of the lowest states in the O"*""*" and 1 -channels along the crossover 
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Figure 5: Connecting the mass spectrum through the crossover, (a) The four lowest 0"*"+ 
states. Empty squares denote the scalar ground state and full squares the glueball. Circles 
correspond to mixed intermediate states, (b) The two lowest 1 states. 



is shown in Fig. ^ The mass of the hghtest state of the theory, the 0"'""'' ground state, 
shows a dip as the crossover is traversed. We have measured this mass on lattices of 
size L = 24, 36 and found no lowering of the mass with increasing lattice size. This 
means that essentially the infinite volume mass is reached and that the physical ground 
state mass definitely stays non-zero at the critical coupling. This identifies the transition 
unambiguously as a crossover It is interesting to note that the mass of the third 

excited state, corresponding to the glueball on the confinement side, does not change 
through the crossover region until well into the Higgs region. This is yet another man- 
ifestation of its decoupling, this time it proves stable under variation of another scalar 
coupling parameter, the hopping parameter Ph- 

In the 1 channel we show only the lightest state and the first excitation, as our 
basis is smaller in this channel. The lowest state clearly becomes lighter as one changes 
from the confinement region to the Higgs region. 

When the parameter is varied through the crossover some states change their 
ordering in mass, so that we felt unable to assign distinct symbols to some of the states 
in Fig. ^(a). In these cases the question arises how one identifies the individual states. 
A useful criterion is to look at the composition of the mass eigenstates in terms of the 
original operators used. As an example, we show the variation of the operator content 
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Figure 6: Operator content through the crossover in the 0~^~^ -channel, a) The ground 
state, b) The glueball. Circles denote the maximal an. from R/L contributions, squares 
those from C and triangles represent the admixture of P^^\ 

of the 0"*""^ ground state and glueball through the crossover region in Fig. ^. The figure 
exhibits a remarkable phenomenon concerning the mixing of operators. As one moves 
from the confinement region towards the Higgs region, the scalar ground state picks up 
more contributions from purely "gluonic" operators. The mixing of scalar and gluonic 
operators is strongest in the middle of the crossover, at the critical coupling (3'^. This is 
not unexpected given the fact that here the fiuctuations are strongest. Although there is 
now a significant gauge contribution to the lowest state, it still remains clearly dominated 
by scalar operators. Moving into the Higgs region the gauge contribution declines again 
and settles at some low value. Looking at the composition of the glueball, however, 
we do not see the analogous feature happening. In contrast, the glueball maintains its 
almost pure gauge composition throughout the crossover beyond the critical coupling, 
with equally low overlap with the Polyakov loop. Only as one moves clearly into the Higgs 
phase does the glueball give up its purely gluonic nature, consistent with the decoupling 
observed from the behaviour of its mass. 

Given the intrinsic interest of this decoupling phenomenon it is worth investigating 
it in more detail. To do so we take the lowest 13 states in the 0^"*" sector and display for 
each state the maximum overlap onto scalar {R or L) operators as well as the maximum 
overlap onto the Wilson loop (C) operators. We do so separately for our extreme confining 
and Higgs coupling values as well as for three values of (3h that correspond to traversing 



26 



1 

Max(R,L) 
0.5 



1 

Max(R,L) 
0.5 



■ XB 

_^ x ■ X 

J > 



1 


_l ' 1 


' L 


Max(R,L) 


x 




0.5 


f 






■ X ^ 

X X 

X ■ ■ 

1 , 1 


X " 





r 



1 


J ' 1 ' L 


Max(R,L) 


■ 




■ 


0.5 


■ 




■ X 




X - 





X 


1,1,1 






1 


J ' 1 ' L 


Max(R,L) 






X 


0.5 


X X 








X 

" ''>S< X X" 





■ 

"1 , 1 , r 




0.5 1 



1 ^ 
It 



_j L_ 



Confinement 
|8h=0.3526 



T 



jgH=0.3531 



_l ' 1 


L 




* 










i8H=0.3538 


i 

1 1 


# _ 

1 





j8n=0.3541 



1 ^ 



# 



j8n=0.3554 
Higgs 







0-5 1 Max(C) 



Figure 7: Variation of the operator content of the 0"*"+ (left) and 2++ (right) spectrum 
through the crossover as described in the teoct. 
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the crossover region from near the edge (on the confining side) to near the pseudo-critical 
value, j3'^. The corresponding plots are shown on the left of Fig. ^ (Errors are not shown 
since they are small.) Those points shown as crosses correspond to states with substantial 
overlap onto fiux loop operators. If the states had random scalar-gauge content, the points 
would be scattered over the whole plot. On the other hand, if there was a total scalar- 
gauge decoupling, then one or the other coordinate of each point would be zero. It is clear 
from Fig. that the latter is close to being the case in the confining phase and that this 
persists, although more weakly, as we traverse the crossover region. In the Higgs phase, 
however, any decoupling effect is clearly much weaker. Note the increasing number of 
states with a substantial flux loop overlap, as we traverse the crossover towards the Higgs 
region. This is a manifestation of the increasing instability of the confining fiux tube that 
we shall discuss in more detail below. 

Also in Fig. 0, on the right, we show similar plots for the 2++. In this case we show 
the errors because several of them are large. The reason for such large errors is that we 
have two almost degenerate states whose order is inverted in some measurements, relative 
to other measurements. That is to say, any point with large errors is definitely not an 
approximate state of the system and should hence be disregarded. If we do disregard 
those "states" then we see that we have a picture very similar to the one we obtained for 
the 0^"*": a marked scalar-gauge decoupling in the confining region that only begins to 
weaken at the very centre of the crossover region. Indeed there seems to be some remnant 
of it even in the Higgs phase. 

6 Flux-tube decay 

When comparing the spectrum on the Higgs and confinement side of the first-order 
phase transition at small scalar coupling we discussed the different interpretations that 
the mass governing the exponential fall-off of the Polyakov line operator has: on the 
Higgs side it can be identified with the volume-independent mass of a two-l^ scattering 
state while on the confinement side it rises linearly with the lattice size and represents a 
fiux loop. These different behaviours represent the different dynamics in the two regimes. 
In this section we try to answer the question how the physics changes from the former 
interpretation to the latter. We know that the two regimes are analytically connected, 
and so are the mass eigenstates of the Higgs and confinement regimes. Hence, there must 
also be a smooth transition in the projection of the Polyakov line operator. 

First, we note that we lose the signal for the mass aMp as we move into the crossover 
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region from either side. At the end-points in both regimes we do see plateaux in the ef- 
fective mass, but they are lost when moving closer towards /3'^. The situation is not 
improved by increasing the statistics. Second, let us recall that the Polyakov loop oper- 
ator has an increasing projection onto the scalar ground state towards the middle of the 
crossover. These two observations may be given a physical interpretation. On the confine- 
ment side of the crossover we saw already a small non-vanishing expectation value for the 
Polyakov loop operator. Nevertheless, the mass extracted from the exponential fall-off of 
its correlator still exhibits the linear dependence on the lattice size and hence represents 
a flux loop winding around the lattice, despite the fact that our model permits screening, 
due to the presence of matter fields. As we move to the critical point the Polyakov loop 
develops a large expectation value which signals the loss of linear confinement, and the 
fact that flux loops no longer exist. The question then is, what happens to the flux loop 
as the parameters are changed? The increasing overlap of the Polyakov loop with the 
Q++ ground state suggests a natural explanation: As we move out of the confining region 
the flux loop becomes more and more unstable, as it is increasingly easy to pair-produce 
scalars that screen the flux. The 0^^ ground state is one possible decay product, just 
based on quantum numbers. A stronger mixing between this state and the flux loop may 
thus indicate the increasing instability of the latter. 

In order to elaborate on this point we consider the overlap of the Polyakov loop with 
all the 0"'""'' eigenstates. If the suggested picture is true then all these states are possible 
decay products. Of course, there are other decay products as well, like pairs of 1 
particles etc. For a first qualitative test of our picture it is however sufficient to consider 
just the single particle states. Let us now treat the flux loop like a decaying resonance. 
As long as it is stable, its two-point correlation function has a pole at its mass. When the 
resonance becomes unstable the pole is shifted away from the real axis. The real part of 
the pole is then equal to the sum of the masses of the decay products, weighted by the 
branching ratios for the particular decay channel. In analogy to this consideration we 
now define an "effective flux loop mass" by the sum of the masses of the 0"*"^ eigenstates, 
weighted by their maximal overlap onto a Polyakov line operator. 



In the same manner we construct the average square of the flux loop energy, (Ep). From 
these we then obtain the "width" of the flux loop as the root mean square fluctuation 



These two quantities are plotted through the crossover region in Fig. ||. Indeed, the 
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Figure 8: Effective energy and decay width of the flux loop for the decay into single 
particle O"*""*" states. 

decay width defined in this way is nearly zero on the confinement side and then increases 
smoothly as it approaches the critical coupling. It then maintains its higher value on the 
Higgs side indicating the instability and non-existence of flux loops on this side of the 
transition. 

We conclude that, although the spectra are very similar on both sides of the critical 
for large scalar coupling, we have clear evidence for a pronounced change in the 
dynamics as one moves from the confining to the Higgs region. In particular it is possible 
to demonstrate the growing instability of flux loops through the transition. We remark, 
however, that the concept of a flux loop only exists in a finite volume, and thus this 
analysis does not have an infinite- volume limit. Nevertheless we would expect the same 
phenomena to happen if one were to investigate the physical flux tube between two static 
external charges at a given separation and then change the parameters from the conflning 
regime to the Higgs regime. 

7 Summary and conclusions 

In this paper we have extended our previous work on the (2 + 1) -dimensional SU(2) 
gauge-fundamental Higgs model in various ways. Firstly, we have added the 2'^'^ mass 
spectrum to our previous calculations of the 0''""'" and 1 spectra. Secondly we have 
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added operators based on Polyakov loops. If we have confinement then these loops will 
project onto chromoelectric flux loops ("torelons") that close through one of the periodic 
spatial boundaries. In that case the lightest mass extracted from the correlation function 
of such an operator will be proportional to the lattice spatial size, and its mass density 
will give the confining string tension (up to corrections, the leading one of which is 
universal and known). This allows us to probe directly the confining properties of the 
theory at various points in the phase diagram. At the same time we have extended our 
basis of operators in the mass spectrum calculations; in part by the inclusion of these 
torelon operators and also products of them. This provides us with a much better control 
of finite-volume effects. We have in addition performed the calculations not only for the 
small scalar self-coupling at which our previous simulations were carried out, but also for 
a large value of the scalar coupling where no phase transition separates the symmetric 
and Higgs phases. Here we have performed calculations through the crossover region that 
separates these two phases. In all these calculations we have used only a single value of 
the inverse gauge coupling, but one for which our earlier calculations assure us that 
lattice corrections to the continuum limit are already small. 

As well as extending our previous calculations, we have also shifted the emphasis of 
our physics interest away from the Elcctroweak Theory at finite temperature to dynam- 
ical phenomena that may be relevant to that other theory with gauge fields and charges 
in the fundamental representation: Quantum Chromo dynamics. 

So one major focus of our interest has been confinement. At small scalar self-coupling 
we find that the torelon mass increases (approximately) linearly with its length and hence 
that the symmetric phase is confining. Of course we know that, just as in QCD, flux tubes 
will be broken through the pair creation of fundamental charges, but within our errors 
we see no sign of this. The string tension is almost identical to its pure gauge value. 
By contrast, in the Higgs phase, we see what perturbation theory tells us to expect: the 
Polyakov loop operator couples to a 1^ 1^-scattering state instead of to the now non- 
existent flux tube. Indeed, apart from the Higgs scalar, W VF-scattering states are all 
we see in the 0+"^ and 2++ spectra. More intriguing is the situation for large scalar self- 
coupling. Here we once again have a symmetric phase that shows linear confinement and 
a Higgs phase where linear confinement is absent. However, now there is no first-order 
transition between the two phases, only a smooth crossover. This raises the question 
of how the dynamics can smoothly interpolate between the two phases. We provided 
evidence for the following picture. In the pure gauge theory a flux loop joined upon 
itself through a boundary is a stable state (of the finite volume Hamiltonian) . With 



31 



fundamental charges present such a state can decay. Far from the crossover such decays 
are strongly suppressed and the flux loop will be almost stable: a resonance with a very 
narrow decay width. As we enter the crossover region the decay width becomes rapidly 
larger: so the pole in the complex energy plane moves rapidly away from the real axis. 
As (initially) a secondary effect, the real part will also change, because of the newly 
enhanced intermediate states, and so the mass will shift as well. At some point, near 
the centre of the crossover, the pole is so far from the real axis that one completely 
loses a particle interpretation. At this point it is no longer useful to talk of confining 
flux tubes. We have found it convenient to use a periodic flux loop on a finite volume as 
a probe of confinement, but one can clearly use a finite fiux line between fundamental 
sources/charges equally well. 

Our second major interest has been the approximate decoupling between scalar and 
gauge degrees of freedom which we observed in our previous work. There wc found that 
in the symmetric phase there are states which arc almost entirely gluonic and the masses 
of these states coincide with the glueball masses of the pure gauge theory. Wc now find 
that this occurs not only in the sector, but also in the 2++ sector, and at large scalar 
self-coupling as well as at small. Moreover this phenomenon only gradually disappears as 
wc move through the crossover to the Higgs phase. Indeed, there is some evidence that a 
remnant of this phenomenon survives even in the Higgs phase. Amongst other things this 
tells us that the mixing between glueballs and states composed of fundamental charges 
is strongly suppressed, at least as long as we have some remnant of linear confinement. 
It also tells us that this is a very robust phenomenon. Is it robust enough to extend 
to QCD? There it would support the idea that glueball-quarkonium mixing should be 
suppressed and that we should search experimentally for glueballs at the masses found 
in the corresponding pure gauge theory. These are important questions but we need 
some understanding of the underlying dynamics if we are to move beyond the realm of 
conjecture. 

In addition to the above two topics, we also clarified the unusually strong finite- 
volume effects that we saw in our previous work. It turns out that these are operator 
artefacts. In fact, the true bound states show weak finite volume effects. However, if one 
uses smeared operators it is difficult to avoid components that are essentially (smeared) 
Polyakov loops. These components introduce the torelon into the measured spectrum 
and, of course, the mass of this varies rapidly with volume. This is a generic problem 
when one has a theory with fundamental charges, and the simplest way to deal with it 
is to include smeared torelon operators in the operator basis. After diagonalisation one 
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can explicitly determine which states have dominant torelon components. 

Finally we point to the dense spectrum of states we have found in the Higgs phase 
at large scalar coupling. This is in a phase without confinement, and presumably the 
dynamics is not being driven by the gauge coupling becoming large at large distances. 
It is therefore tempting to conclude that what we have here is a generic example of a 
spectrum produced by a strong scalar self-coupling. 

We conclude with a brief comment as to how these calculations need to be improved. 
The most obvious lacuna concerns multiparticle scattering states. We have not attempted 
to identify these, except in the very obvious case of the Higgs phase at small scalar self- 
coupling. There is some evidence that our operators, which are naively designed to project 
only onto single particle states, usually do so. However, this really needs to be checked 
explicitly, especially in the strong-coupling Higgs phase where, if it turns out not to be 
the case, it could completely alter our interpretation of the spectrum. The inclusion of 
scattering states is also necessary if we are to perform a proper analysis of flux tube 
decay through the crossover region. Technically the problem is simple. For example, 
suppose t) is an operator that has a very good projection onto the lightest scalar 
with momentum p (which will be one of the discrete set of allowed momenta in our finite 
volume). Then the operator 02 (^) = <l>{+P,'t)<Pi~P,'t) (vacuum subtracted if necessary) 
will be a good operator for the two-scalar scattering state of zero total momentum and 
finite relative momentum. At least this will be so for volumes that are not too small. 
In similar ways one can build other trial multiparticle operators and add them to the 
operator basis. However, to be confident that one is correctly isolating multiparticle 
states one needs to do an accurate finite volume analysis that shows the continuum 
cut gradually forming as one increases the volume. This would be a substantially larger 
calculation than the one in this paper. It would also be very interesting to further pursue 
the scalar-glue decoupling; in particular to see if it might be related to the OZI rule in 
QCD. 
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